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ABSTRACT 

We have used multi-epoch long-baseline radio interferometry to determine the proper motion and 
orbital elements of Algol and UX Arietis, two radio-bright, close binary stellar systems with distant 
tertiary component s. For Algol, we re fine the proper motion and outer orbit solutions, confirming 
the recent result of Zavala et al. (20101 that the inner orbit is retrograde. The radio centroid closely 
tracks the motion of the K1V secondary. In addition, the radio morphology varies from double-lobed 
at low flux level to crescent-shaped during active periods. These results are most easily interpreted 
as synchrotron emission from a large, co-rotating meridional loop centered on the K-star. If this is 
correct, it provides a radio-optical frame tie candidate with an uncertainty ±0.5 mas. For UX Arietis, 
we find a outer orbit solution that accounts for previous VLBI observations of an acceleration term 
in the proper motion fit. The outer orbit solution is also consistent with previously published radial 
velocity curves and speckle observations of a third body. The derived tertiary mass, 0.75 solar masses, 
is consistent with the Kl main-sequence star detected spectroscopically. The inner orbit solution 
favors radio emission from the active K0IV primary only. The radio morphology, consisting of a 
single, partially resolved emission region, may be associated with the persistent polar spot observed 
using Doppler imaging. 

Subject headings: astrometry, stellar dynamics, binariesxlose, radio continuum:stars 



1. INTRODUCTION 

Multi-epoch phase-referenced radio interferometry is 
a powerful method to determine astrometric properties 
of radio- loud objects in a variety of astrophysical set- 
tings. Very long baseline interferometry (VLBI), with 
baselines exceeding 10,000 km, can provide astrometri c 



accuracies between 10 and 100 /ias (Pradel et al. 2006), 
allowing direct distance measurements to many gaiac- 
tic and extragalactic radio sources with un precedented 



2009 


Deller et al. 


2009 


Dzib et al. 


2010) 



g., Reid et al. 



metric ~VLBl a pplicatio ns" inc lude tests of general rela- 



tivity ( Fomalont et al.||2009|), prec i se determination o f 
Earth's mertial frame ( |Petrov||2007| |Petrov et aTf 2009 ) , 
and sear ches for extra-sola r planets in late-type stellar 
systems (Bower et al. 2009). 

In this paper, we report on astrometric VLBI studies 
of two well-studied active stellar systems, Algol and UX 
Arietis. Both systems contain a very active K sub-giant 
tidally locked in a short period binary system with a dis- 
tant third companion. In addition, both are relatively 
close (< 50 pc), so the angular scale of the inner orbit is 
several times the resolving beam of a global VLBI array. 
This allows tracking the motion of the active component 
within the inner binary orbit. By imaging the radio-loud 



component over many epochs, it is possible to determine 
all astrometric parameters (parallax, proper motion, in- 
ner and outer orbital elements) with high accuracy. In 
addition, the resulting maps provide a direct probe of the 
geometry of the coronal structure. Finally, if the solution 
is sufficiently accurate, the position of active star can be 
inferred on the radio image, allowing a tie between the 
optical re ference frame and the radio ICRF with sub -mas 
accuracy (Lestrade et al. 1999 Bourda et al.| 2010). 

Algol (HD 19356, B8V+K21V, 29 pc), the prototype 
of the eponymous Algol binary class, has been well- 
studied since its identification as an eclipsing binar y sys- 



tem more than two centuries ago ( Goodricke|l~783 ) . The 
cooler subgiant secondary fills its Roche lobe, causing 
episodic mass transfer onto a n accretion disk surroundin g 
the main sequence primary ( |Richards fc Albright |p93| . 
The subgiant is chromospherically and coronally a ctive 



with frequ ent s trong flares obse rved at radio (e.g., Mu- 
I tel et al.| [19981 jRichards et al.| |2Q Q3j), ultraviolet (e^T 
Istern et al. 19 95h, and X-ray "Te.g., Ottmann & Schmitt 
|l996||!Schmitt &; Eavata|1999 ) wavelengths. The close bi- 
nary pair has an orbital period of 2.83 days and is tidally 
locked. A distant tertiary companion (F1V, P = 1.86 yr) 
is in an eccentric orbit oriented nearly perpendicular to 
the inner binary orbit. 
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Although most of Algol's orbital parameters for both 
inner and outer or bits were well-determined spectroscop- 
ically decades ago (iHill et al. 11971 Bachmann & Hershey 
T975] |Stein & Beardsley||l977|), 



several parameters were 
either undetermined or ambiguous until more recent di- 
rect imaging became possible. The inner binary's ori 



entation was first determined by Lestrade et al. ( 1993 1 
who observed the active K-star's radio emission at both 



quadratures using a global VLBI array. Pan et al. ( 1993 ), 
using optical interferometric data, determined that the 
orientation of the outer orbit was nearly perpendicular 
to the inner orbit, confirming the position angle inferred 
from speckle data ( jBonneau 1979). 

Determination 01 the inclination of the inner orbital 
plane has proved more problematic. Csizmadia et al. 
( |2009[ ) imaged Algol using both an optical interferome- 
ter and VLBI. They found that the inner orbit was pro- 
grade, i.e. that the line-of-sight component of the orbital 
ang ular momentum poi nted toward the observer. How- 
ever Zavala et al. ( 2010 ) imaged both the inner and outer 
binaries at multiple epochs using the NPOI 6-element in- 
terferometer and found that the inner o r bit is retrograde. 



In addition, they found that Pan et al. (1993 1 's determi 
nation of the outer orbit ascending node longitude was 
in error by 180°. 

Interpreting the physical properties of a stellar radio 
source critically depends on accurately registering the 
radio morphology with its location within the stellar sys- 
tem. Algo l was the first st ellar system ever imaged using 
VLBI ( |Clark et al~1|1975[ ). They found that during an 
exceptionally large radio flare (S ~600 mJy), the source 
characteristic size was comparable to the overall angu- 
lar size of the inner binary system. The inferred high 
brightness temperature and broad spectrum were con- 
sistent with gyrosynchrotron emission from mildly rela- 
tivistic electro ns in a coronal magnetic field. Subsequent 
VLBI images flMutel et al.||1985| |Lestrade et aT1|1988[ 
Mutel et al.||1998[ |Massi et al.||2002| |Csizmadia et a ' 



guiar resolution and astrometric accuracy of these ob- 
servations taken separately prevented an accurate regis- 
tration with th e stella r components. However, recently 



Peterson et al. (2010) observed Algol at six epochs us- 
ing a high-sensitivity global VLBI array at a higher fre- 
quency (15 GHz) than past observations. These data 
allowed much better sensitivity and angular resolution. 
They found that the radio structure consisted of a large 
coronal loop that was centered on the active K star, 
aligned along the inner binary axis, and co-rotating. In 
this paper we combine this dataset with archival phased- 
referenced VLBI observations made since 1983. This 
combined dataset allows a global solution for all astro- 
metric parameters. 

UX Arietis (HD 21242, G5V+K0IV, 50 pc) is an active 
double-lined spectroscopic binary system. The primary 
K subgiant shows activity at r adio through UV wave- 
lengths very similar to Algol ([Carlos fc Popper 1971 
Trigilio et al.|1998[|Mutel et al.|1998||Torriceln-Ciamponi 



mi 
19; 



eFaIT l998; Lang & Willson 1988: Bu ccino fc Mauas|2009 
Ekmckci 2010). There is also chromospheric activity as- 
sociated with the secondary, possibly caused by episodic 
mass-transfer associated with Roche-lobe overflow from 
the primary ( Hu enemoerder et al.|[l989l |Gu et aLl|2002 ; 
,s & 



Hatzes (1991) used Doppler imaging to detect a large, 
stable polar spot on the K primary, as well as tran- 
sient spots at intermediate latitudes, which app ear to be 
preferential ly oriented toward the G companion ( Rosario 



et al. 2008). As with Algo l, VLBI imaging o f the ra- 



dio corona of UX Arieties ( Mute l et al. 1984 Beasley 



fe Gudel||2000[ |Massi et al. 2002 ) shows structure com 
parable to the inner binary separation and a brightness 
temperature consistent with gyrosynchrotron emission. 

A third component has long been suspected in the 
UX Arietis system, based on both radial velocity anoma - 
lies (|Duemmler fe Aaruml|2001| |Glazunova et al.||2008|), 
speckle observations of a distant third body (e.g., |jVlcAl 
ister et al.|[T987] |Hartkopf et~H1|2000"l |Balega et aT 



2006 ) , and apparent acc eleratio n in proper motion stud- 



les 



2006 



2009) confirmed this basic picture, but the limited an- jLestrade et al 



.estrade et al.||1999| |Boboltz et al.||2003[ |Fey et al 
). In this paper, we combine new multi-epoch VLBI 
measurements with archival VLBI data to solve for all 
astrometric components including a tertiary component 
in a long-period orbit. 

2. OBSERVATIONS 

Algol and UX Arietis were observed with a global VLBI 
array at six and four epochs respectively between June 
2008 and October 2009. The array comprised ten 25 m 
telescopes of the Very Long Baseline Array (VLBA|^[ 
the 100m Green Bank Telescope (GBT), and the 100m 
Effelsberg telescope (EB). Observations were performed 
at 15.4 GHz in dual-polarization mode with a bandwidth 
of 128 MHz in each polarization. The array synthesized 
beamsize (~ 0.5 mas ) was considerably smaller than the 
projected angular separation of the inner binary orbit of 
either system, so we were able to map motion within the 
inner binary. In addition, we used both previously pub- 
lished and archival VLBA observations to supplement 
our data, allowing highly accurate global astrometric so- 
lutions over a time internal spanning more than 25 years. 
We used the 'nodding' phase-referencing te c hnique 



|lont et al.p roS 



1990| |Beasley fc Conway|[l995| |Foma 
switching rapidly between a compact 
extragalactic iC RF so urce with well-determined coordi- 
nates (Ma et al. 1998) and the target star. For the Al- 
gol observations, we used two phase calibrators: The 
angularly close primary phase calibrator ICRF source 
J031301.9+41200 (A6» = 1.0°) and a secondary calibra- 
tor, J031049.8+381453 (A0 = 2.8°). We alternated be- 
tween two-minute scans on Algol and one-minute scans 
on the primary phase calibrator, with additional one- 
minute scans of the secondary calibrator every hour. For 
the UX Arietis observations, we alternated between scans 
of UX Arietis and ICRF source J033630.1+321829 (a.k.a. 
NRAO 140, A9 = 4.2°) using a 90-second cycling time. 

In addition to these observations, we included two 
other radio astrometric measurements for these stars. 
First, we included sever al previously published radio 



interferometer positions (jLestrade et "ah 1993 1999 



Boboltz et al.||2003| |Fey et al.||2006| ). SeconTJ, a num- 
ber of phase-referenced observations of both Algol and 



UX Arietis exist in the NRAO VLBA archive ( Sjouwer 
man et aL]|2004 ). We calibrated and imaged these data 



Aarum Ulvas & Engvold 2003a; Ekmckci 20l0p. |Vogt"fc| 



1 The National Radio Astronomy Observatory is operated by 
Associated Universities Inc., under cooperative agreement with the 
National Science Foundation. 
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and present the resulting positions in this paper for the 
first time. All VLBI data analyzed for this paper are 
summarized in Table [T] 

3. ANALYSIS 

3.1. Data Calibration and Imaging 

VLBI visibilities from all epochs were calibr ated and 
image d using the NRAO AIPS software package (Greisen 
2003), except for the 1983.5-1994.4 UX Arietis positions 
tha t were analysed using the SPRINT software package 



as radial velocity data (Ducmm ler &: Aarum 2001 Hill 



et al. 


1971 


Hill et al. 


1 


terential positions from 



93| IGudel et al. 1999) and dTF- 



(cf. |Lestrade et al.||1999[ ). We followed standard VLBI 
amplitude and delay/ rate calibration procedures. Geo- 
metric delays introduced by small errors in the predicted 
values of the Earth Orientation Parameters (EOP) were 
corrected by downloading tables of the measured EOP 
for the time of each observation. We also corrected iono- 
spheric and tropospheric delays — these are discussed in 
section 13.31 below. 

The phase calibrator sources used in the obser vations 
are compact extragalactic core-jet sources (e.g. Mutel 
et al. (1998)). We made self-calibrated images of the 
calibrators (SNR ^5000:1) which we used to recompute 
the fringe solution. In principle, the core positi on of a 
calibrator c ould shift with frequency. However, Foma- 



lont et al. (2011) observed four ICRF calibrators from 
8-43 GHz and found the core positions are coincident to 
within 0.02 mas. Hence, we have assumed that the core 
centroids of the calibrators have a negligible shift from 
8-15 GHz. 

For the Algol observations, we also corrected for source 
motion within the observing period. Since the radio 
emission in Algol is associated with the K sub-giant 
(Lestrade et al. 1993), it may move nearly 2 mas dur- 
ing a 12-hour observation, especially if the observation is 
centered near one of the eclipses. Images made without 
correcting for this motion would result in a smeared-out 
source along the direction of motion. Using the already 
well-determined orbital elements of the inner binary in 
Algol, we introduced a time-varying phase-tracking cor- 
rection to the visibilities (using AIPS task CLCOR) 
equal to the K star's position offset in the inner binary 
orbit. 

In order to determine radio centroid positions, we made 
Stokes I images using target visibilities that had been 
corrected from linearly interpolated delay-rate and phase 
solutions of calibrator scans bracketing each target scan. 
We then self-calibrated these images in order to improve 
image fidelity and better discern the source structure. 
If the source was unresolved, we fit a Gaussian to the 
centroid of the emission (AIPS task JMFIT) and used the 
Gaussian FWHM/2 as the position uncertainty estimate. 
For epochs with a resolved double structure we found 
the centroid of each radio lobe and used the midpoint 
between the lobes as the radio position for that epoch. 
In both cases, the position uncertainties were increased 
if the calculated position uncert aint y due to tropospheric 
and ionospheric delays (section 3.3) exceeded the image- 
based estimate. 

3.2. Parameter Fitting 

We performed a global least-squares minimization of 
the orbital and astrometric parameters designated 'vari- 
able' in Tables [3] and [I] using our VLBI data as well 



We began with a coarse grid search, varying each pa- 
ramcter by small steps over a reasonable range of val- 
ues and finding the minimum x 2 -value i n the resulting 
multi-dimensional parameter space. We used the Nelder 
Mead simplex (a.k.a. 'amoeba') algorithm (Press et al. 



1992[ ) to refine the best-fit parameter values from the 
coarse grid search. This minimization scheme finds the 
direction of steepest descent without needing to calcu- 
late partial derivatives, which is especially useful when 
solving transcendental orbital equations. It also uses an 
adaptive stepsize in order to converge to within a very 
small tolerance of the minimum chi-square solution in 
the parameter space. 

3.3. Astrometric accuracy 

An ideal interferometer with project ed baseline length 
B ha s a positional accuracy given by (Thompson et al. 
19861 

° 2tt SNR B ' { ' 

where SNR is the signal-to-noise ratio of the target 
source and A the observing wavelength. This ideal es- 
timate, of order several /Ltas for centimeter-wavelength 
VLBI, is almost never realized in practice since instru- 
mental and intervening atmospheric effects degrade the 
phase of the incoming signals. 

Phase-referenced astrometric observations using mod- 
ern VLBI arrays typically have very small antenna posi- 
tion and timing uncertainties. Neither of these uncer- 
tainties contribute significantly to the astrometric ac- 
curacy so long as the switching time between calibra- 
tor and source is smaller than the coherence time. The 
two most important sources of position uncertainty result 
from path lengt h fluctuations in the troposphere and the 
ionosphere (e .g^ Ftamalont et al.| 1995[ Ros k, Re"id||2005 
Pradel et al.| 2006 ). A good phase-referencing scheme 
(i.e. with a snort enough cycle time) should compensate 
for these delays, but for sub-mas astrometry they can 
still be a source of error and should be addressed. 

The daytime ionospheric delay is dispersive (r cx A 2 ) 
and is highly dependent on time of day. At 10 GHz, 
the daytime zenith ionospheric delay is quite large, typi- 
cally 1 ns, corresponding to an excess path length 30 cm. 
Fortunately, since 1998 there have been several global 
databases of ionospheric delay maps available for online 
download. These maps are generated every two hours 
and a re derived from observed GPS delays (Ros et al. 
2000[ ). We have corrected all po st-1998 VLBI da ta tor 
Ihese delays (AIPS task TECOR). |Ros et al.|(|200io| ) found 
that the residual ionospheric uncertainty after this cor- 
rection is less than ±0.15 ns (4.5 cm) at 8.4 GHz on 
global baselines. 

The tropospheric delay is non-dispersive, and consists 
of a dry component well predicted by the local atmo- 
spheric pressure, and a highly variable wet component 
whose value can exceed 10 cm excess path length. Both 
GPS delay maps and water vapor radiometers have been 
used to measure and correct for the wet component (e.g. 
Snajdrova et al.|2006[ ) , reducing the uncertainty to about 
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TABLE 1 
VLBI Observing Log 



Exp. Code Nr. epochs Dates Freq. (GHz) Cal. Source Array Ref. 



BM044 
BM074 
BM109 
BM267 



Algol 

1 1995.309 8A 0313+412 VLBA + Y27 

3 1997.303 - 1997.308 8.4 0313+412 VLBA +Y27 + EB 

3 1999.042 - 1999.058 8.4 0313+412 VLBA + Y27 + EB 

6 2008.265 - 2008.628 15.4 0313+412 VLBA + GBT + EB 



IMutel et al. 



1 1998 I 



Peterson et al 



] (|2010} 



UX Arietis 





9 


1983.568 - 


1994.405 


8.4 


0326+277 


(a) 


BB032 


1 




1994.833 


8.4 


0326+277 


VLBA 


BB049 


6 


1995.875 - 


1995.885 


8.4 


0326+277 


VLBA 


BG097 


1 




2001.128 


8.4 


0326+277 


VLBA + EB 


BM140 


4 


2001.728 - 


2001.739 


8.4 


0326+277 


VLBA + EB 


BP157 


1 


2009.638 - 


2009.797 


15.4 


NRAO140 


VLBA + GBT + EB 



Lestrade et alJ7l999 1 
ae 



Beasley k, OvAli\ l |2UU0| 



|Massi et al] |2002[ ) 
This paper 



a VLBI arrays varied with epoch, see |Lestrade et a.l.| ( |199Q] ) 

3 cm. However, although all VLBA sites measure the 
local atmospheric pressure and hence the dry compo- 
nent accurately, they are not equipped with water va- 
por radiometers, so there is no on-site measurement of 
wet component delay. Rather, the correlator applies a 
seasonally-based average value, which can be in error by 
several cm on a given day, especially for humid sites such 
as St. Croix. 

For phase-referenced observations on a single baseline, 
these path delay un certainties result in a position uncer- 
tainty given by (cf. |Reid et aL]|1999[ ) 



c St 
A 



tan(Z) sec(Z) ■ SZ ■ 9 , 



(2) 



where Z is the mean zenith angle, SZ is the zenith angle 
difference between calibrator and target, 5tq is the zenith 
path length uncertainty, and is the baseline angular 
resolution. 

In Figure [l] we show a representative plot of the posi- 
tional uncertainty expected from ionospheric and tropo- 
spheric delay uncertainties as a function of frequency for 
a global VLBI array using the nodding technique at mean 
zenith angle 70° and target-calibrator angular separations 
l°and 4°. We have also assumed a tropospheric uncer- 
tainty ±0.1 ns (3 cm), and mean daytime ionospheric 
uncertainty of ±0.1 ns (3 cm). Nighttime ionospheric de- 
lay fluctuations are much smaller (~10% daytime values) 
and are unimportant at frequencies above a few GHz. 

It is clear that at 15 GHz, the tropospheric contribu- 
tion dominates, resulting in position uncertainties 0.12 
mas at 1° calibrator-target separation and 0.5 mas at 4° 
separation. At 8.4 GHz, the daytime ionospheric and the 
tropospheric contributions are nearly equal, resulting in 
a total uncertainty 0.18 mas at 1° separation and 0.68 
mas uncertainty at 4° separation. Unless given other- 
wise, we use these uncertainties as formal uncertainties 
for all VLBI measurements in Table [l] 



4. RESULTS 

Table [2] lists radio centroid positions for both sources 
determined at all VLBI observing epochs. The global 
least-squares astrometric solutions computed from these 
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Fig. 1. — Position uncertainty vs. frequency for baseline length 
10 4 km, 70°zenith angle, tropospheric delay uncertainty ±0.1 ns 
(3 cm) and daytime ionospheric delay uncertainty ±0.1 ns (3 cm). 
The solid and dashed lines show position uncertainties for target- 
calibrator separations of l°and 4°, respectively. Position uncer- 
tainties due to tropospheric delay uncertainty are shown by the 
green lines, the daytime ionosphere contribution is shown by red 
lines, and the nighttime ionosphere contribution is blue. The 
15.4 GHz position uncertainties are dominated by the troposphere, 
while at 8.4 GHz, the daytime ionospheric contribution is compa- 
rable to the tropospheric contribution. 

positions are listed in Tables [3] and |4j The correspond- 
ing differences between observed positions at each epoch 
and those calculated using the global solution are given 
in Table [5] where the column labeled a is the number 
of standard deviations between the observed and model 
positions. For both solutions, the agreement is excellent, 
with only a few epochs marginally above 2cr. 



4.1. Algol 
4.1.1. Outer orbit 

Figure [3] shows the position of the radio-loud secondary 
(Algol B) on the sky, after correction for proper motion 
and parallax. The model trajectory is shown as gray 
line, while the observed and predicted positions at each 
observing epoch are shown as blue errorbars and red x 's, 
resp ectively. Our outer orbital solution agrees very well 
with|Zavala et al. (2010), including the orientation of the 
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TABLE 2 
Observed radio centroid positions 



JD Mid-observation <j> a S a <5 a a erg Freq. Ref. 

2,400,000+ UT Inner geocentric geocentric heliocentric heliocentric mas mas GHz 



Algol 



47629 


375 


1989/04/12 


21 


00 





28 


03 


08 


10 


127004 


40 


57 


20 


34098 


03 


08 


10 


125682 


40 


57 


20 


32477 





75 


1 


50 


8 


4 


1 


47632 


302 


1989/04/15 


19 


14 





30 


03 


08 


10 


127162 


40 


57 


20 


34201 


03 


08 


10 


125973 


40 


57 


20 


32589 





75 


1 


50 


8 


4 


1 


47633 


448 


1989/04/16 


22 


15 





70 


03 


08 


10 


126933 


40 


57 


20 


33808 


03 


08 


10 


125797 


40 


57 


20 


32200 





75 


1 


50 


8 


1 


1 


47636 


333 


1989/04/19 


19 


59 





71 


03 


08 


10 


127074 


40 


57 


20 


33958 


03 


08 


10 


126073 


40 


57 


20 


32364 





75 


1 


50 


8 


4 


1 


49831 


333 


1995/04/23 


19 


59 





23 


03 


08 


10 


129554 


40 


57 


20 


33252 


03 


08 


10 


128722 


40 


57 


20 


31680 





65 





85 


8 


1 


2 


50559 


354 


1997/04/20 


20 


30 





13 


03 


08 


10 


130422 


40 


57 


20 


32512 


03 


08 


10 


129472 


40 


57 


20 


30925 


1 


15 


1 


55 


8 


1 


2 


50560 


354 


1997/04/21 


20 


30 





48 


03 


08 


10 


130524 


40 


57 


20 


32603 


03 


08 


10 


129621 


40 


57 


20 


31021 





85 





75 


8 


4 


2 


50561 


354 


1997/04/22 


20 


30 





83 


03 


08 


10 


130358 


40 


57 


20 


32413 


03 


08 


10 


129504 


40 


57 


20 


30838 


1 


00 


1 


15 


8 


1 


2 


51194 


583 


1999/01/16 


01 


59 





67 


03 


08 


10 


128550 


40 


57 


20 


34164 


03 


08 


10 


125901 


40 


57 


20 


33939 


1 


65 


1 


15 


8 


1 


2 


51197 


583 


1999/01/19 


01 


59 





72 


03 


08 


10 


128552 


40 


57 


20 


34124 


03 


08 


10 


125843 


40 


57 


20 


33815 





60 





75 


8 


4 


2 


51200 


583 


1999/01/22 


01 


59 





76 


03 


08 


10 


128433 


40 


57 


20 


33862 


03 


08 


10 


125672 


40 


57 


20 


33469 


1 


10 


1 


30 


8 


1 


2 


51889 


271 


2000/12/10 


18 


29 





95 


03 


08 


10 


130700 


40 


57 


20 


34500 


03 


08 


10 


129269 


40 


57 


20 


35259 


1 


00 


7 


00 


8 


1 


3 


52212 


000 


2001/10/29 


12 


00 





50 


03 


08 


10 


134300 


40 


57 


20 


33500 


03 


08 


10 


134977 


40 


57 


20 


35027 


2 


70 


3 


10 


8 


4 


4 


54084 


375 


2006/12/14 


21 


00 





50 


03 


08 


10 


133525 


40 


57 


20 


32636 


03 


08 


10 


131934 


40 


57 


20 


33302 


3 


00 


3 


00 


5 





5 


54563 
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UX Arietis 
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Note. — (a) Heliocentric positions were calculated using the parallax given in Tables [3] and H] (b) Coordinates were measured with 
respect to primary phase calibrators : 0313+412 (•03:13:01.962129, + 41:20:01.18353'), 0326+277 (03:29^7.669425, +27:56:15.49901), NRAO140 
(03:36:30.107609, +32:18:29.34226) ||Ma et al.|1998||Fey et al. | 2004|). (c) Previously published star coordinates have been correct ed for different 
assum ed c alibrator position s if needed. Id) Keterences : 1. [Lestrade et al.|j 1993ft, 2. Positions from VLBA archival data, 3. |Boboltz et al. 
J2003|), 4 . JFev et al.) ( |2006| t , 5. [Csizmadia et al.|f2009| , 6. This paper, 7. |Lestrade et al.J jl984| |1999[ l and unpublished positions, 8. |Beasley 
& Uudcl (2000), position calculated trom V LldA archival data 9. [Franciosini et al. ( 1999|, positions calculated from VLBA archival data. 
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TABLE 3 
Algol orbital element solutions 



Parameter 


Symbol 


Value 


unit 


Param.type 


Parallax 


n 


34.7 ±0.6 


mas 


Variable 


Primary mass 


m A 


3.70 


M 


Fixed 


Secondary mass 


m B 


0.79 


M e 


Fixed 


Tertiary mass 


mc 


1.51 ±0.02 


M e 


Calculated 


R.A. proper motion 


Ma cos 5 


2.70 ± 0.07 


mas yr 1 


Variable 


Dec. proper motion 


US 


-0.80 ± 0.09 


mas yr -1 


Variable 


Fiducial" R.A. 


ao 


03:08:10.132409 ± 0.7 


mas 


Variable 


Fiducial a Declination 


o 


40:57:20.3353 ± 0.6 


mas 


Variable 


Fiducial epoch 


JD 


2452212.02 


JD 


Fixed 






2001.82553 


year 




Inner binary 


Period 


Pi 


2.867329 


days 


Fixed 


Eccentricity 


ei 


0.000 




Fixed 


Inclination 


h 


99 ± 17 


deg 


Variable 


Longitude of ascending node 


Hi 


48 ± 20 


deg 


Variable 


Lonigtude of periastron 


0J1 


270 


deg 


Fixed 


Semi-major axis 


Ol 


2.26 


mas 


Fixed 


Time periastron 


Ti 


2441773.49 


JD 


Fixed 






1973.24638 


year 




Outer binary 


Period 


P-2 


679.5 ± 0.3 


days 


Variable 






1.86 


years 




Eccentricity 




0.16 ± 0.02 




Variable 


Inclination 


i'2 


85.5 ± 1.4 


deg 


Variable 


Longitude of ascending node 


n 2 


130.7 ± 3.5 


deg 


Variable 


Longitude of periastron 




312.0 ± 1.4 


deg 


Variable 


Semi-major axis 


0-2 


95.4 ± 0.5 


mas 


Variable 


Time periastron 


T 2 


2446930.0 ± 3.2 


JD 


Variable 




1987.36551 


year 





Center of mass of triple system at fiducial epoch 



7 



TABLE 4 

UX Arietis orbital element solutions 



Parameter 


Symbol 


Value 


unit 


Par am. type 


Parallax 


n 


19.90 


mas 


Fixed 


Primary mass 


m A 


1.10 


M 


Fixed 


Secondary mass 


m B 


0.95 


M Q 


Fixed 


Tertiary mass 


m c 


0.75±0.01 


Mq 


Calculated 


R.A. proper motion 


f^a cos 8 


44.96 ± 0.13 


mas yr 


Variable 


Dec. proper motion 




-102.33 ±0.09 


mas yr — 1 


Variable 


FiduciaP R.A. 


ao 


03:26:35.38386 ± 1.2 


mas 


Variable 


FiduciaP Declination 


So 


28:42:54.2755 ± 0.8 


mas 


Variable 


Fiducial epoch 


JD 


2451544.5 


JD 


Fixed 






2000.0 


year 




Inner binary 


Period 


Pi 


6.4378553 


days 


Fixed 


Eccentricity 


ei 


0.000 


Fixed 


Inclination 


h 


59.2 


deg 


Fixed 


Longitude of ascending node 


Qi 


82 ± 34 


deg 


Variable 


Longitude of periastron 


Wl 


180.0 


deg 


Fixed 


Semi-major axis 


ai 


1.71 


mas 


Fixed 


Time periastron 


Ti 


2450642.00075 


JD 


Fixed 






1997.52705 


year 




Outer binary 


Period 


P2 


40548.6 ± 70.2 


days 


Variable 






111.02 


years 




Eccentricity 


e2 


0.77 ± 0.01 




Variable 


Inclination 


»2 


93.3 ± 0.6 


deg 


Variable 


Longitude of ascending node 


02 


58.9 ± 0.5 


deg 


Variable 


Longitude of periastron 


W2 


274.9 ± 0.8 


deg 


Variable 


Semi-major axis 


«2 


648.0 ± 0.7 


mas 


Variable 


Time periastron 


T-2 


2451664.9 ± 34.3 


JD 


Variable 






2000.32964 


year 





a Center of mass of triple system at fiducial epoch 

Algol C orbit, which differed from earlier determinations 
(see Zavala et al. for discussion). We note that the de- 
rived mass of the tertiary component (f .57 ± 0.01 Mq), 
although determined to high accuracy, depends on the 
masses of A and B, which we have assumed are exactly 
known. 

We can use the outer orbit solution to compute the or- 
bit of Algol C relative to the AB center of mass and com- 
pare with optical observations. The model AB-C orbit is 
shown in Figure pi green line) along with optical obser- 
vations of Algol (J tabulated in the Fourth Catalog of In- 
terferometric Measure ments of Binary Stars (FCIMBS) 
( |Hartkopf et al. 2010). The observed positions are in 
good agreement with the model orbit and predicted po- 
sitions at each observing epoch. 

4.1.2. Inner orbit 

Figure [4] shows the orbit of Algol B in the inner bi- 
nary center of mass frame, along with the mean observed 
(blue errorbars) and model (red x's) positions at the 
six observing epochs observed at 15 GHz (experiment 
BM267). We have not displayed the 8 GHz positions, 
since their position uncertainty is twice as large as the 
15 GHz points and provide little additional constraint of 
the inner binary orbit. The orientation of the inner bi- 
nary, £1 = 48° ±20°, a grees with the earlier VLBI result of 
Lestrade et al.l (I1993I) (52° ± 5°) and the more recent op- 



tical determination on |Zavala et al.| ( |2010[ ) (47.4° ±5.2°). 
We also find i = 99° ± 17", conhrming the conclusion of 
Zavala et al. (20101, who found that the inner orbit is 




Ac* cos(i>) (mas) 



Fig. 2. — Algol C (tertiary) orbit with respect to the AB center 
of mass determined from VLBI global astrometric solution (green 
line) along with speckle interferometer observations listed in the 
Fourth Catalog of In terferometric Measurements of Binary Stars 
| |Hartkopf et al.|2010} (blue errorbars) and model positions for the 
corresponding dates (red x's). The green arrow is positioned with 
its tail at the point of periastron, and shows the direction of orbit. 
The green circle is at the ascending nodal point. 



retrograde (i.e. i > 90°, contrary to the previous deter- 
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TABLE 5 

VLBI — MODEL POSITION DIFFERENCES 



Aa cos(J) (mas) 



Fig. 3. — Orbital path of the Algol inner binary center of mass 
(light gray line), with VLBI observations (blue errorbars) and 
model positions of the radio-loud star (Algol B, red x's). The 
model orbit of Algol B is shown at the time of each observation 
(dark gray ovals) for reference. The gray arrow shows the orbital 
direction for the AB center of mass, and the gray circle is the as- 
cending node. 




Fig. 4. — Algol inner binary system, showing model orbit (gray 
line), and observed (blue errorbars) and predicted (red x's) radio 
centroids at six epochs between 2008.27 — 2008.63 . The gray 
arrow shows the direction of orbital travel. The ascending node of 
the orbit lies toward the upper-left in the figure. 



urination of Csizmadia et al. ( 2009 ). 



4.1.3. 



Algol's radio morphology: Co-rotating coronal loop 
centered on KIV star 

Algol's radio morphology consists of a double-lobed 
structure oriented normal to the inner orbital plane dur- 
ing quiescent states, an d a crescent or loop st ructure dur- 
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optical maps cannot be made directly, since the optical 



ing active periods (see Peterson et al.| 2010 and supple- 
mentary materials). The registration between radio and 



positio n uncertainty ellipse (5x8 mas, Perryman et al. 
( |1997| ) is larger than the inner binary orbit. However 
our 15.4 GHz observations show that the radio centroid 
mirrors the predicted motion of the KIV secondary over 
six epochs to an uncertainty ±0.5 mas. 

It is possible that the radio centroid is systematically 
displaced from the K star and is following its trajectory 
precisely, but this seems unphysical. Occam's razor leads 
us to conclude that the correct registration is that the 
radio centroid is coincident with the K star center. In 
that case, the lobe structure is straddling the active KIV 
secondary, since the lobe separation (~1.0 mas) is the 
same as t he K star's angular dia meter (1.1 mas). As dis- 
cussed by Peterson et al. (2010 and online supplement), 
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the radio morphology and registration is consistent with 
a co-rotating, plasma-filled coronal loop structure emit- 
ting synchrotron radiation. The loop is oriented in the 
direction of the primary, which may imply magn etic in- 
teraction w ith the primary's accretion disk (e.g., Retter 
et~aLl[2005 ) 



X-ray observations (Chung et al. 2004) show that the 
centroid of the quiescent emission orbits with a semi- 
major axis 15% smaller than the K star — i.e. it is offset 
toward the center of mass of the inner binary. The radio 
positions in FigureHlappear likewise shifted inward of the 
K star orbit. Our global astrometric solution found that 
the semi-major axis of the centroid of the radio emission 
is shifted inward by 0.2 ± 0.8 mas. It is worth noting that 
the shift in the radio centroid is flux-dependent, moving 
inward toward the top of the co-rotating loop when the 
high-flux events fill the loop. 

4.2. UX Arietis 

Unlike Algol, the orbital elements of a third component 
in the UX Arietis system were not previously known, al- 
though the presence of a third component is not unex- 
pecte d — companions i n clos e binaries are quite com- 



Tokovinin et al. (2006) found that 63% of spec- 



troscopic binaries [SBJ that they surveyed had at least 
a tertiary companion. For SB's with a period less than 
3 days, the fraction with companions is 96%, suggesting 
that the shorter period systems exchanged angular mo- 
mentum with their companions, shortening their orbital 
periods. 

In the UX Arietis system, a third component had been 
inferred by detectio n of non-linear proper motion (e.g., 
Boboltz et al. | 2003[) and spectral lines of a po ssible third 
component (Aarum Ulvas & Engvold 2003a). However, 



two quite different models: a 10.7 yr circular orbit and a 
21.5 yr highly eccentric orbit (Duemmler & Aarum 2001 ). 
Therefore, it was necessary to undertake a comprehen- 
sive parameter search, constrained by th e VLBI posi 
tions, third c omponent radial velocity data (|Duemmler fc 
Aarum|2001||Aarum Ulvas &: En gvold 2003a; Glaz unova 



the putative third component's orbit was poorly deter 
mined. The only published orbital solution consisted of 



et al. 2008p , a nd optical inte rfcromctric observations in 
the FCIMBS ( |Hartkopf et al | |2010[ ). The search was 
complicated by the relatively high proper motion of UX 
Arietis (/i ~100 mas yr" 1 ), which magnifies a small frac- 
tional proper motion uncertainty over decades into a 
large aggregate position shift. 

4.2.1. Outer orbit 

Figure [5] shows the radio position of UX Arietis on the 
sky plane, from epoch 1980.0 (origin) to 2009.8. The blue 
circles are observed positions relative to 1985.0. The un- 
dulating shape of t he trajectory, previously interpreted 
as an accelerati on (|Lestrade et al. 1999 Boboltz et al. 
2003 Fey et al. 2006)) results from the reflex motion ot 
the inner binary m a 111-yr period outer orbit. The po- 
sition of the radio component in the frame of the outer 
binary center of mass in shown in Figure [6| The gray 
line is the predicted position of the inner binary center 
of mass, with blue errorbars indicating observed posi- 
tions and red x 's the model positions after correction 
for proper motion and parallax. 

We next compare the model solution with FCIMBS in- 
terferometer observations. Figure [7] shows the model po- 
sition of the tertiary component with respect to the inner 
binary center of mass (green line). The red and corre- 
sponding blue symbols are the model and observed posi- 
tions of the tertiary component as listed in the FCIMBS. 
Note that for all FCIMBS observations prior to 2002, we 
have applied a 180°correction to the catalogued position 
angle. Finally, Figure [8] shows the predicted radial veloc- 
ity of the inner binary center of mass (gray line) and ter- 
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Fig. 5. — UX Arietis radio position projected on the sky plane 
since 1980.0. The undulating trajectory results from proper motion 
combined with reflex motion in the AB-C outer binary. 
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Fig. 6. — Position of UX Arietis primary relative to triple system 
center of mass, showing model (gray line), and predicted (red x's) 
and observed (blue +'s) positions. The arrow shows the orbital 
direction and is positioned with its tail at the point of periastron. 
The gray circle shows the ascending nodal point. 
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Fig. 7. — UX Arietis C (tertiary) orbit with respect to the AB 
center of mass determined from VLBI global orbit solution (green 
line), along with speckle interferometer observations listed in the 
Fourth Catalog of In terferometric Measurements of Binary Stars 
| |Hartkopf et al,|2010|l (blue errorbars) and model positions for the 
corresponding dates (red x's). The green arrow is positioned with 
its tail at the point of periastron, and shows the direction of orbit. 
The green circle is at the ascending nodal point. 



radial velocity observations flDuemmler & Aarum 


2001 


Massarotti et al. 


2008; 


Glazunova et al. 2008|. in gen- 



although the outlier at epoch 1991 25 from the HIPPAR- 
COS catalog ( |Perryman et al. 19971 is not consistent with 
the orbit solution. 
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Fig. 8. — UX Arietis outer orbit radial velocities jDuemmler ~Ez\ 



|Aarum|200l"l|Massarotti et al.|2d08| jGlazunova et aI.[2U08| | 

curves based on the outer orbit solution are sfiown tor tl 



Model 
~t!he inner 

binary center of mass (gray line) and tertiary (green line), along 
with observed velocities for the tertiary (green squares) and AB 
center of mass (blue circles). 



4.2.2. Inner orbit 

Since both inner binary comp onents of UX Arietis ex- 
hibit c hromospheric activity (Aarum Ulvas & Engvold 
2003a), it is unclear whether the radio emission originates 
from a single component, both components or perhaps 
an accretion region between the two stars. The orbital 
elements are well-characterized by spectroscopic observa- 
tions ex cept for the orientation (Q ) and the inclination 
(i ^60°, [Duemmler fc Aarum|2001 1, which is degenerate 
about 90 reflection with respect to radial velocity curves. 
Hence we fixed the other orbital parameters and, as with 
Algol, used only the 15 GHz observations to solve for for 
Q and i. 

We confirm that the KOIV primary is the source of the 
radio emission in UX Arietis by testing the quality of the 
fit assuming the alternatives: that the emission is either 
from the secondary, or from a location between the two 
components. The alternative assumptions both yielded 
unacceptable fits (xt = ?-6 an d xt = 2-8, respectively), 
while the solution for radio emission from the primary 
was very good, yielding a good fit at the 99.7% confidence 
level. 

Figure [9] shows the best-fit orbital solution overlaid 
with contour maps of the radio emission at each epoch. 
We obtain a tentative solution for Q (82° ±34°). Our 
data cannot constrain the solution in i, as can be seen by 
the fact that the timing of our observations were all very 
close to the nodal points in the inner binary orbit. Fix- 
ing the inclination at 59.2° (Duemmler & Aarum 20011 
produces a slightly better fit than the K. V .-degenerate 
value (i — 120.8°), thus we use the former value in our 
final solution. 



4.2.3. Source structure 

For all four VLBI observations made at 15 GHz, 
the self-calibrated radio images are amorphous, with 
an overall size comparable to the restoring beam 
(~0.4x0.8 mas). This contrasts with earlier published 



images of UX Arietis flMutel et al. 1985 Beasley & 



Gudel||2006| |Ros fc Massi||20U7p ~wmch~lhow structure 



on an angular scale comparable to the binary separa- 
tion (1.7 mas). However, the previous observations were 
at lower frequencies (5 GHz and 8.4 GHz) so it is pos- 
sible that the interferometer phases were corrupted by 
uncorrected ionospheric delay fluctuations, especially for 
observations made prior to 1998, when GPS-based iono- 
spheric corrections were unavailable (cf. section 3.3). Al- 
ternatively, there may be frequency-dependent extended 
structure with a steep spectral index that is not de- 
tectable at higher frequencies. 

5. DISCUSSION 

5.1. Mass and spectral type of UX Arietis tertiary 
component 

The UX Arietis outer orbit solution can be used to de- 
rive the mass of the tertiary component, m c = 0.75 ± 
0.01 Mq, where the formal uncertainty is contingent on 
fixed values for the inner binary masses, as given in Ta- 
ble |4| This mas s corresponds to a spectral type Kl main 



sequen ce star (Zombeck 1990| . Duemmler & Aarum 
( |2001| estimated the mass of UX Arietis' tertiary com- 
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Fig. 9. — UX Arietis inner orbit solution overlaid on radio contour maps at epochs (a) Aug 21, (b) Sep 15, (c) Oct 15, and (d) Oct 18, 
2009. Orbits and predicted positions of the KIV primary (blue line, circle) and secondary (red line, circle) are shown. The blue arrow 
shows the orbital direction and position of the ascending node. The contour levels of the radio images are 10%, 30%, 50%, 70% and 90%. 
of peak values 0.7, 5.2, 3.8, and 1.4 mjy beam - in panels a, b, c, and d respectively. The total flux density at each epoch is (a) 0.9 mjy, 
(b) 9.8 mjy, (c) 7.5 mjy and (d) 2.8 mjy. 



ponent between 0.30 < m 3 < 0.46 M Q , for their cir- 
cular and elliptical orbit solutions respectively, although 
they reject the smaller mass, since it would be an M 
star too faint to be re sponsible for the spectral lines as- 
cribed to the tertiary. Aarum Ulvas & Engvold (2003a I 
fit the continuum spectrum of the UX Arietis system 
and found that after correction for the primary and sec- 
ondary spectra, the derived color index is best-fit to a 
K5 main-sequence star, although they cast doubt on this 
identication, since the implied radius (0.83 solar radii) is 
in the range of a Kl-type star. Our derived mass agrees 
with the hotter type, so perhaps starspots on the ter- 
tiary are infl uencing its color index as is susp ected with 
the primary | | Aarum Ulvas fc Engvold||2003b| ). 



5.2. Comparison with polar loop models 



Peterson et al. (2010) recently described a polar loop 
model tor Algol in which the double-lobed radio compo- 
nents are associated with foot points of an active region 
of meridional magnetic field lines originating at the star's 
magnetic poles. A simil ar geometry was already sug- 
gested for this system by Franciosini et al.| ( |1999[ ). We 
now consider whether the single radio component seen 
in the UX Arietis system can be interpreted within this 
model. 

A large, stable polar spot has been detected on the K 
prima ry of UX Arietis ( Vogt fc Hatzesfl99T Elias et al. 
|1995| ). These polar spots are a common phenomenon: 
Ktarspots on magnetically active cool stars preferentially 
appea r near the poles (e.g. |Hatzes et al.fl9 96 Bruls et al. 
1998). This preference for high latitudes may be due to 
rapid rotation, which leads to the Coriolis force dominat- 
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ing over buoyancy forces in the dynamics of magnetic flux 
tubes. As a consequence, flux tubes in the stellar convec- 
tion zone migrate nearly parallel t o the axis of rotation 
and t hus surface at high latitudes (Schuessler & Solanki 
1992|. Hence, magnetically active stars with rapid rota- 



tions to calculate proper motions and orbital elements of 
the inner and outer orbits in both stellar systems 



For Algol, we confirm the early result of Lestradc et al 
( 1993 1 that the radio centroid closely tracks the motion 



tion exhibit magnetic flux eruption at high latitudes and 
polar starspots. 

The VLBI images of UX Arietis indicate that the ra- 
dio centroid lies in the hemisphere facing away from the 
G star. This is consistent with observations of a simi- 
lar late-type binary , HR1099, in which doppler imagery 
(Donati et al. 19921 and X-ray monitoring indicate that 
the active regions lie in th e hemisphere facing aw ay from 



ot the KIV secondary. Furthermore, the radio morphol- 
ogy, which varies from double-lobed at low flux level to 
crescent-shaped during active periods, is consistent with 
synchrotron emission from a large, co-rotating merid- 
ional loop centered on the K-star. If this correct, it pro- 
vides a radio-optical frame tie candidate with a precision 
± 0.5 mas. We also refine the proper motion and outer 
orbit solutio ns, confirming the recent result of |Zavala| 



the inactive companion. Audard et al. (20011 suggest 



that this is because tidal locking of the rotation of the K 
star alters the internal dynamo in such a way that strong 
activity on the hemisphere facing the inactive star is sup- 
pressed. 

Since the inner binary orbital plane is inclined at 30° 
to the observer's line of sight (i ~ 60°), if the magnetic 
axis is close to perpendicular to the inner binary's or- 
bital plane, only one pole will be visible. The polar loop 
model posits that the radio emission arises above both 
polar regions, but in this case, the far side emission will 
be largely occulted by the K star. This geometry may 
explain the absence of two radio centroids. 

6. SUMMARY 

We have analyzed more than 20 years of phase- 
referenced VLBI observations to determine accurate 
(sub-mas) positions of the radio centroids from two ac- 
tive close binaries, Algol and UX Arietis, both of which 
have distant tertiary companions. We used these posi- 



et al. (120101) that the inner orbit is retrograde. 

For UX Arietis, we find an tertiary orbit solution that 
accounts for previous VLBI observations of an acceler- 
ation term in the proper motion fit, as well as radial 
velocity curves and speckle observations. The dynam- 
ical mass, 0.75 ± 0.01 solar masses, supports the iden- 
tification of the tertiary with a Kl main sequence star, 
consistent with third-body color index measurements of 
Aarum Ulvas & Engvold (2003a). The inner orbit solu- 
tion favors emission from the active K primary only, in 
the hemisphere facing away from the G star. The ra- 
dio morphology, consisting of a single, partially resolved 
emission region, may be associated with the persistent 
polar spot observed using Doppler imaging. 



The authors would like to thank Dr. Bob Zavala of the 
U.S. Naval Observatory Flagstaff Station for his helpful 
comments and insights, without which our paper would 
have been much less complete. 
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